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The Yellow Excitonic Series of Cu 2 Revisited by Lyman Spectroscopy 
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We report on the observation of the yellow exciton Lyman series up to the fourth term in 
CU2O by time-resolved mid- infrared spectroscopy. The dependence of the oscillator strength on 
the principal quantum number n can be well reproduced using the hydrogenic model including 
an AC dielectric constant, and precise information on the electronic structure of the Is exciton 
state can be obtained. A Bohr radius ai s =7.9 A and a ls-2p transition dipole moment (J,i s -2p= 
4.2 eA were found. 
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The concept of exciton, a composite quasi-particle in a 
semiconductor consisting of an electron-hole pair bound 
by Coulomb attraction, is well established. 1-3 Based on 
the effective-mass approximation, the internal electronic 
structure of an exciton can be described by a two-particle 
Schrodinger equation leading to a series of bound states 
similar to that of atomic hydrogen. In the case of an ex- 
citon, however, the Rydberg energy is reduced by several 
orders of magnitude because of the screening effect of the 
crystal lattice. It should also be noted that an exciton is 
the electronic elementary excitation of a system of many 
interacting electrons, so that a simple two-body picture 
has fundamental limitations. 4 The many body aspect of 
excitons has recently attracted renewed interest, related 
to issues such as the interpretation of excitonic lumines- 
cence, 5 the origin of excitonic nonlinearities, 6 and the oc- 
currence of a Bose-Einstein condensation (BEC) phase. 

Recent advances in laser technology have made new 
spectroscopic tools available. In particular, mid, or far- 
infrared spectroscopy, often called 'THz spectroscopy,' 7 
allows to probe intraband kinetics of electrons and holes, 
the collective motion of electron hole ensembles, as well 
as the transitions between the different terms of an ex- 
citonic series. With excitonic Lyman spectroscopy, it is 
also possible to probe optically inactive excitons such as 
spin triplets, and invisible excitons lying outside the light 



cone 



8-10 



The yellow excitonic series in CU2O is a textbook ex- 
ample of Wannier excitons. A very clear hydrogenic se- 
ries is observed up to the 10 th term in linear absorption 



spectroscopy. 
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Because of the positive parity of the 



valence (symmetry I\jT) and conduction (symmetry Ig ) 
bands, the yellow series starts with the n=2 term. The 
Is yellow series exciton level is split due to electron-hole 
exchange interaction into a forbidden singlet paraexciton 
(I2 ; 2.021 eV at 2K) and a triply degenerate orthoex- 
citon (Tg" ; 2.033 eV at 2 K). The spin configuration of 
paraexcitons is a purely triplet state and thus optical 
transitions to the ground state are strictly forbidden to 
all orders. Consequently, the Is paraexciton has a very 
long lifetime and it has long been recognized as a promis- 

*E-mail address: gonokami@ap.t.u-tokyo. ac.jp 




Fig. 1. Schematic energy diagram for the yellow series Is and np 
exciton states and their Lyman series transition in CU2O. 



ing candidate for the observation of exciton BEC. 15,16 
The resonance energies of the n >2 terms obtained by 
linear absorption measurement can be well reproduced 
using the hydrogen formula, thus providing the basic pa- 
rameters to characterize hydrogenic wave functions. The 
resonance of the Is exciton state, however, has a signif- 
icant deviation from the position extrapolated from the 
higher Rydberg lines (n >2). Such a deviation is known 
as the central cell correction of the Is exciton state. It is 
caused by the non-parabolicity of the electron and hole 
bands at large wave vector due to the small Bohr radius 
of Is excitons (comparable to the lattice constant 17 ) and 
to dynamical screening. There still remain some ambigu- 
ities on the precise electronic structure of the Is exciton 
even if this information is crucial to determine the sta- 
bility of the exciton BEC phase. 

In this paper, we report on the observation of the ex- 
citon Lyman series corresponding to transitions from a 
population of cold Is orthoexcitons, created by a reso- 
nant two-photon process, to higher p states up to the 5 
term. Figure 1 schematically illustrates the time-resolved 
exciton Lyman absorption process occurring in our ex- 
periment along with an energy diagram of the yellow 
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series excitons in C112O. The split between ortho and 
paraexcitons is 12 meV. At a low density, the orthoex- 
citon relaxes to the lower energy paraexciton state via 
a process involving the emission of transverse acoustic 
phonons, with a relaxation time of a few nanoseconds. 18 
Recent experiments with a continuous wave probe light 
allowed to measure the Lyman signal associated with the 
ls-2p transition of the paraexcitons. 19 ' 20 Time-resolved 
measurements with femtosecond pulses revealed the pres- 
ence of both ortho and paraexcitons. 8,9 

A systematic analysis of the oscillator strength and 
line shape of the Lyman series allows to examine the de- 
tailed features of Is excitons such as the wave function 
of the relative motion of electron-hole pair, the trans- 
lational mass, and the distribution function of Is exci- 
tons. For such an analysis, it is necessary to obtain first 
the Is excitons with a narrow energy distribution. Since 
orthoexcitons can be generated by resonant two-photon 
excitation, 21 we can apply a phase space compression 
scheme with ultrashort laser pulses, in analogy with the 
two-photon excitation of biexcitons in CuCl. 22 With this 
scheme, we can prepare orthoexcitons inside a very small 
phase-space, with an energy spread much smaller than 
the laser pulse bandwidth and the thermal energy of the 
lattice. In a one-beam excitation configuration, the or- 
thoexcitons have an additional common kinetic energy 
of the the order of a few /zeV which is gained from the 
finite photon momentum. 

The details of the experimental setup and procedure 
have been described in previous papers. 8-10 A naturally 
grown single crystal has been cut in the form of a platelet 
shape with a 220 /im thickness along the (100) plane. The 
sample is cooled by contact with a copper block main- 
tained at liquid Helium temperature. A tunable pump 
source with a wavelength centered around 1220 nm and 
a 500 fs pulse width (8 nm spectral width) is used. The 
pulse length is shaped from a 100 fs pulse by using a 
pair of gratings. The excitation pulse with 1 [ii energy 
is focused onto a 3xl0 4 yum 2 area of the sample. The 
excitation light is linearly polarized and its orientation is 
set to maximize two-photon absorption. A mid-infrared 
tunable source with a wavelength centered around 10 ^m 
is used for the weak probe pulse. The transmitted probe 
light is analyzed by a monochromator (0.1 meV spectral 
resolution) followed by a HgCdTe detector. The Cu 2 
crystal shows a suitable transparency window between 
112 and 137 meV, allowing the observation of transitions 
from the Is orthoexciton to the np states. 

Figure 2(a) shows a linear absorption spectrum in the 
spectral region of np orthoexciton states. In Figure 2(b), 
an induced absorption spectrum obtained under reso- 
nant two-photon excitation of Is orthoexcitons at a low- 
density excitation (pump pulse energy 3 mJ/cm 2 ) mea- 
sured after a pump-probe delay of 5 ps is shown. Under 
this excitation condition, a low exciton density was ob- 
tained, with a value estimated to be n exc iton ~10 15 cm~ 3 
(See below). Four lines are observed and are positioned 
at 116, 129, 133, and 136 meV respectively. One can see 
that the spectral features obtained in Fig. 2(b) are in 
very good agreement with the ones shown in Fig. 2(a) 
if the energy shift corresponding to the Is orthoexciton 
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Fig. 2. Linear absorption spectrum (a). Induced absorption spec- 
tra of the Is orthoexciton generated instantaneously by two- 
photon excitation (b), and the Is paraexciton measured at 4 
ns after the excitation (c). 



line -Ei sor ift, o =2033 meV 24 is taken into account. Thus 
we can assign the series of Fig. 2(b) to the transition 
from the Is to the np orthoexcitons (n=2,3,4,5), i.e., the 
Lyman series. Their peak positions, widths, and area of 
absorption are summarized in Table 1. 

It should be noted that the energy split between the 2p 
and 3p orthoexctions is close to that between the Is ortho 
and paraexcitons. The ls-2p paraexciton line should be 
carefully distinguished from the ls-3p orthoexciton line 
as described in our previous paper. 9 Figure 2(c) shows a 
transient absorption spectrum measured at a 4 ns pump- 
probe delay. The ls-2p orthoexciton signal is weaker and 
the higher energy narrow lines are not present due to the 
decay of orthoexcitons into paraexcitons, while a rela- 
tively broad line appears around 128 meV. This new line 
is assigned to the ls-2p transition of the paraexcitons. 

The lincwidth of the transitions of the exciton Lyman 
series has a systematic n-dependence. A similar behav- 
ior has been observed in np absorption spectra and has 
been discussed in terms of an n-dependent LO phonon 



Table I. Peak position, peak widths, and normalized integrated 
surfaces of the ls-np (rt=2,3,4,5) lines. AEi 3 — np is the energy 
estimated from the Is and np absorption lines. 
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Fig. 3. Lincwidth of ls-np lines (o) and rap absorption lines (A). 

scattering rate from np to Is state. 25 Figure 3 shows the 
linewidths of the ls-np transitions (o). The total width 
of the lines has several contributions: the energy distri- 
bution of Is excitons, the intrinsic linewidths of the Is 
state, and the width of the final np states. The finite 
energy distribution of the Is excitons, f\ s , contributes 
to the broadening because of the difference in the mass 
parameters between the Is and np excitons schemati- 
cally shown in Fig. 1. The widths of the Lyman lines 
can be expressed as: (AEi s _ np ) 2 ~(AE np ) 2 +((mi s /m2 P - 
l)As) 2 - The intrinsic linewidth of the Is exciton state is 
much narrower (7i s <C0.1 meV) than the spectral reso- 
lution of the present experiment and is neglected in our 
analysis. We estimate the linewidth of the np states by 
fitting the linear absorption data lines with asymmet- 
ric Lorentzian curves following the model of Toyozawa. 26 
The obtained result is also plotted in Fig. 3 (A). From 
the difference of the width between the ls-4p and the 4p 
lines (<0.1 meV), the upper limit of the energy spread 
of initial Is orthoexcitons, f\ s , is estimated to be 0.36 
meV/(3m e /1.7m e -l)<~0.5 meV, which is narrower than 
the energy spread expected for the 4.2 K thermal equi- 
librium exciton gas, 0.7 meV. This indicates that the 
initial super cooled Is orthoexcitons are still at a tem- 
perature below that of the lattice even at a delay of 5 
ps. 

The wave functions of np states are well described by 
a hydrogenic model, since the spectral positions of the 
np exciton states (n >2) are well described with a Ry- 
dberg series £"„=17525-786/n 2 cm -1 . 13 However, the Is 
exciton state shows a significant deviation from this Ryd- 
berg series which has been discussed by several authors. 
Most recently, Artoni ct al. analyzed the AC Stark ef- 
fect under resonant pumping of the ls-2p transition of 
orthoexciton and used a value /ii s _2 P =6.3 cA (30 D) 
in direct analogy with hydrogen atoms. 27 Kavoulakis et 
al. examined the shrinkage of the Bohr radius of the Is 
exciton state caused by the central cell correction due 
to the non-parabolicity of the valence and conduction 
bands. 17 They took into account the quadratic terms 
of the wave vector dependence of the bands under the 
constraint that the known value of translational mass of 
Is orthoexciton, 3m e must be reproduced. This analysis 
gave a Bohr radius a\ s =5.3 A and a binding energy 
E=164.9 meV. They found that the wave function of 
the Is exciton is well approximated with such a modi- 
fied Is function. 2 Recently, it was pointed out that such 



Fig. 4. (a) Integrated spectral area of ls-3p lines (o) and ls-4p 
lines (A) as a function of area of the ls-2p lines, (b) Integrated 
area of ls-np induced absorption lines normalized by that of ls- 
2p line. +: Experimental results, o: Best fit (4.2 eA). 

a small effective Bohr radius of the Is excitons reduces 
the overlap between the Is and 2p exciton wave func- 
tions, yielding a revised dipole moment of 1.64 eA. 23 
The systematic analysis of the exciton Lyman absorp- 
tion provides a unique opportunity to determine the "un- 
known" Is exciton wave function from a set of projections 
to the known np states through dipole matrix element 

f^ls—np G <^ Ui s \z\u n p ^> . 

The integrated area of the exciton Lyman absorption 
line ls-np can be expressed as follows: 

L-n P ^)dE = ^lEl^tl^l (!) 

where Ei s ^ npi ii\ s ~ np and e are the ls-np transition en- 
ergy, dipole moment and background dielectric constant 
at the probe frequency respectively. We determine the 
background dielectric constant e ~6 by measuring the 
interference fringe pattern of the absorption spectrum 
in the mid-infrared region. Knowing the density of the 
Is exciton and from the measured integrated absorp- 
tion Si s - n p = J ls _ np Aa(E)dE, we can determine the 
dipole moment /ii s - np . In practice, however, it is difficult 
to estimate the absolute density of the Is exciton from 
the absorbed photon number, which strictly depends on 
the detailed pumping condition especially in the case of 
a two-photon excitation with short pulses. Therefore in 
order to determine the wave function of the Is exciton 
state, we rely on the relative intensity of the integrated 
absorption of the ls-np lines Si s - np / S\ S -2 P , which are 
density independent and insensitive to the pump condi- 
tion. Indeed, in the low-density excitation limit, Eq.(l) 
yields the ratio, S ls - np / S ls - 2p =| Mis-n P /Mis-2 P | 2 

Figure 4(a) shows the integrated spectral area of 
the ls-np induced absorption (n=3,4) as a function 
of the area of the ls-2p transition. By varying the 
excitation intensity, we find that the integrated areas 
of the ls-3p and ls-4p transitions are proportional to 
the integrated area of the ls-2p line, i.e. the Is ex- 
citon density. From the fitting of the result shown in 
Fig. 4(a), we obtain the ratio S , i s _ 3 p/5 , i s _ 2 p=0.26±0.01 
and S'i s _4p/S'i s _ 2 p=0.11±0.01. Figure 4 (b) shows 
the Si s _„p/Si s _2p integrated area of the ls-np lines 
(n=3,4,5) normalized by that of the ls-2p line. The 
crosses show the experimental results. With our anal- 
ysis, we obtain S ls - np / Si s - 2p = |Mis-rap/ ' His-2 P \ 2 x 
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(E ls - np /E ls - 2p ) where 

and Ei s _ n p are the transition dipolc moment and res- 
onance energy of ls-np transition respectively. From the 
observed value of S\ s - np / S\ s ^2p, we obtain a best fit 
value of the Is Bohr radius of ai,5=7.9 A, which gives 
a dipole moment of ^i s _ 2p =4.2 eA. It should be noted 
that the Is Bohr radius extracted from the positions 
of the np exciton series (n>2) is oi s =ll.l A. Accord- 
ing to Eq. (1), the spectrally integrated absorption area 
Si s -2p ~1 meVcm -1 corresponds to an exciton density 
ni s ~6xl0 13 cm -3 under our relatively low excitation 
conditions. 

Equation (1) implies that we can determine the value 
of the transition dipolc moment from the absorption 
strength if we know the density of excitons. To cross 
check the validity of our procedure to determine the 
dipole moment and Bohr radius based on the ratio 
2p , we also measured the induced Lyman ab- 
sorption under one photon excitation at a 600 nm wave- 
length, corresponding to the phonon side band of the Is 
orthoexciton, where we can expect a quantum efficiency 
close to 100 % for photon-exciton conversion. With an 
absorbed photon density of 2xl0 15 photons/cm 3 , we 
measured an induced absorption signal at a 50 ps delay 
S ~53 cm meV in the spectral region from 112 meV 
to 137 meV, which includes the contributions from the 
ls-np orthoexcitons and the ls-2p paraexcitons. At such 
short delay time, we can neglect the decay of the orthoex- 
citons. We obtain /Ui s _2p of 4.8 eA. The value is slightly 
overestimated since the contribution of the ls-3p transi- 
tion is neglected. Thus the obtained value is consistent 
with the estimation from the S\ s - np /S\ s -2 P ratio. 

With a hydrogenic model, the main cause for the devi- 
ation of the Is exciton level from the Rydberg series is the 
dynamic effect on the electron-hole screening. Here, we 
introduce the AC dielectric constant at a frequency close 
to the exciton binding energy. The measured dielectric 
constant is e ~6. The effective Bohr radius is determined 
from the binding energy Eu n and dielectric constant e 
with <xb — 27re 2 /i?bi„eeo. Eu n ^150. 8 meV and e ^6.0 
in turn give a value as ~7.8 A, which agrees with our 
estimate. This indicates that the hydrogenic model with 
only a correction due to screening by the AC dielectric 
constant is suitable to explain the general feature of Is 
excitons in Cu 2 0. 

In summary, we observed the Lyman series up to the 
5p term of the yellow series exciton in Cu 2 using fem- 
tosecond time-resolved mid-infrared spectroscopy. From 
the measured narrow lincwidths of the higher Rydberg 
series, we inferred that the Is orthoexcitons have an ef- 
fective temperature lower than the lattice temperature 
even 5 ps after their creation by direct two photon exci- 
tation. This confirms that phase space compression of the 
resonant two-photon excitation with femtosecond pulses 
is an effective method to create a super cooled exciton 
gas. From a systematic analysis of the Lyman series ab- 
sorption strength, we obtain a Bohr radius of the Is ex- 
citon ai s =7.9 A and a ls-2p transition dipole moment 



Mis-2 P =4.2 cA. The obtained dipole moment enables us 
to determine the exciton density directly from the spec- 
trally integrated absorption area; ni s /5i s _2p=6x 10 13 
cm _3 /meVcm~ 1 . This novel method for the quantita- 
tive measurements of exciton density is important for 
the study of high-density effects in excitons including ex- 
citon Bosc-Einstcin condensation, Auger recombination 
and exciton Mott transition. 
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